Compressibility, zero sound, and effective mass of a fermionic dipolar gas at finite 

temperature 
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The compressibility, zero sound dispersion, and effective mass of a gas of fermionic dipolar 
molecules is calculated at finite temperature for one-, two-, and three-dimensional uniform sys- 
tems, and in a multilayer quasi-two-dimensional system. The compressibility is nonmonotonic in 
the reduced temperature, T/Tf, exhibiting a maximum at finite temperature. This effect might 
be visible in a quasi-low-dimensional experiment, providing a clear signature of the onset of many- 
body quantum degeneracy effects. The collective mode dispersion and effective mass show similar 
nontrivial temperature and density dependence. In a quasi-low-dimensional system, the zero sound 
mode may propagate at experimentally attainable temperatures. 



I. INTRODUCTION 

The study of ultracold atomic systems has received 
much attention in recent years, motivated largely by 
the prospect of realizing novel strongly correlated many- 
body physics. A series of remarkable experimental break- 
throughs have produced an extremely well-controlled 
physics playground [H-Q. The latest breakthrough is 
the ability to associate atoms via a Feshbach resonance 
to form ultracold molecules in the ro-vibrational ground 
state [1, [1], and the JILA group has achieved a ne arly 
degenerate gas of ground state polar molecules [1, Ist-fTTf. 
This is a true milestone in the program, since strong dipo- 
lar interactions between degenerate molecules in the pres- 
ence of an external electric field allow for the design of 
exotic Hamiltonians [H, [l^l and are expected to give rise 
to exciting phenomena including roton softening [14l - ll6l | , 
supersolidity [T7l - [2l| , artificial photons [13] , bilayer quan- 
tum phase transitions f23l| , and multi- layer self-assembled 
chains (23 | for bosonic molecules, and spontaneous in- 
terlayer superfluidity [251 . itinerant ferro electricity p6| . 
Fermi liquid anisotropy |27H29| , fractional quantum Hall 
effect 130J, Wigner crystallization ISll. biaxial nematic 
phase |33| . topological superfluidity j34|, and a Z2 topo- 
logical phase [35| for fermionic molecules. 

Most, if not all, of these novel quantum phases will 
require temperatures on the order of O.lTp or less (with 
Tp the Fermi temperature), which will require further 
experimental advance. Efforts to overcome the current 
obstacles of collisional instability and insufficient cool- 
ing are already underway, though, and given the rate of 
progress in reaching the current state-of-the-art, the fu- 
ture is bright 38]. However, it is worthwhile to ask 
what interesting effects one might be able to observe in 
the immediate future with the temperature on the order 

of Ti;-. 

From a condensed matter perspective, the very idea of 
having a 1 /r^ interaction rather than the usual Coulomb 
interaction is intriguing. Such a system has no paral- 
lel in solid state materials, and study of the previously 
unmotivated problem of many-body physics in a system 
interacting via a potential is still in its infancy. Of 



course, the actual form of the inter-molecular interac- 
tion potential is quite complicated, particularly at short 
range, and the scattering and stability of the molecules 
are sensitive to these details j39| . Further investigation 
is required to allow full quantitative comparison between 
theory and experiment. For our purposes, though, it 
is sufficient to take a interaction with some unde- 
termined short-range cutoff; the low-energy many-body 
physics of a stable, dilute gas should not depend qualita- 
tively on the short-range details. It is then appropriate 
to consider what quantum many-body effects might be 
calculated for a gas with f /r^ interactions that could be 
observed in cold polar molecule experiments. 

In this paper, we calculate the compressibility (or 
equivalently, the ordinary sound dispersion), zero sound 
dispersion, and effective mass at finite temperatures com- 
parable to those achieved in current experiments. These 
quantities should be readily accessible to experiment, and 
we present the first calculations to include thermal ef- 
fects, as well as some trap effects. As these calcula- 
tions are carried out within a leading order perturbation 
approximation in the dipolar coupling constant, we are 
working within the standard weak-coupling theory in the 
sense of the Landau Fermi liquid theory. This is the first 
theory one must do before one does anything else for the 
dipolar systems. 

For dipolar interactions, in contrast to the familiar 
Coulomb case, the weak-coupling regime corresponds to 
the /ow-density limit. This is fortuitous, since current 
experiments cannot achieve high densities without signif- 
icant loss rates. For all practical purposes, the leading- 
order perturbative results are exact at typical experi- 
mental densities. However, lower densities correspond 
to higher T/Tp, necessitating the calculations be done 
for finite temperature. As the density is increased, the 
leading-order perturbative results should remain qualita- 
tively correct as long there is no phase transition to break 
adiabaticity. For very high densities one will enter the in- 
teresting strong-coupling regime and the weak-coupling 
theory will fail in a systematic way, but this regime is 
completely inaccessible at present in dipolar molecular 
systems. 
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In two-dimensional semiconductor-based electron sys- 
tems and in graphene, studying compressibility experi- 
mentally [iol - li^ and theoretically [43, has been an 
important tool for studying quantum many-body effects 
in Coulomb systems. In these two-dimensional condensed 
matter systems with the 1/r Coulomb interaction, the 
Hartree-Fock approximation works remarkably well for 
understanding compressibility, even in the strongly inter- 
acting regime. This is a general result of the frequency 
independence of the compressibility, and not specific to 
the Coulomb interaction. 

We allow for different trap geometries by considering 
uniform three-dimensional (3D), two-dimensional (2D), 
and one-dimensional (ID) g well as a nonuniform 

3D gas in a periodic potential along z. We find that under 
certain conditions the compressibility varies nonmono- 
tonically with temperature. A closely related quantity 
with the same temperature dependence, dEp /dfi (with 
Ep the Fermi energy and the chemical potential), also 
varies nonmonotonically with density, even for tempera- 
tures on the order of Tp. The zero sound speed and effec- 
tive mass also exhibit nontrivial dependencies. However, 
in 3D, propagation of zero sound requires lower temper- 
atures than currently feasible. In 2D and ID, though, we 
find that zero sound propagates undamped even at exper- 
imentally realistic temperatures, assuming that the in- 
termolecular potential behaves like at short enough 
distances. In the low-dimensional situations, the geom- 
etry with dipoles aligned perpendicular to their separa- 
tions is particularly relevant as this is the most stable 
configuration conducive to our approximation of a 
interaction, so we will focus on this case. 

The layout of the paper is as follows: In section HIl we 
review the form of the dipolar interaction in momentum 
space. In section lllll we present analytic low-temperature 
expansions for the self-energy, compressibility, and effec- 
tive mass of a uniform system, as well as numerical re- 
sults for arbitrary temperature. We also present numer- 
ical results for the zero sound mode. In section IIVI we 
numerically obtain the compressibility and low-lying col- 
lective modes for the spatially inhomogeneous case of a 
ID periodic external potential dividing the 3D cloud into 
multiple quasi-2D layers. In section|Vl we discuss experi- 
mental observation of quantum many-body effects in the 
compressibility, and we conclude in section I VI I In Ap- 
pendix A, we briefly compare with the case of Coulomb 
interactions. 



tion r is given by 
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where 9r is the angle between r and E. We will pa- 
rameterize the interaction using the dimensionless ratio 
A = d"^ /r^Ep, where is the average interparticle dis- 
tance. For the different dimensionalities we have 
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where kp^ is the noninteracting Fermi wavevector. For 
fermionic KRb [d — 0.57 Dcbye, m ~ 127 amu) at densi- 
ties around lO^^^^-a [gj, |8|, |g[ (or IO^cto"^ or lO^cm,-^ 
in lower dimensions), A ~ 0.1. 

Strictly speaking, the Fourier transform of V^d (r) has 
ultraviolet and infrared divergences. However, recall that 
anyways this potential is not valid for arbitrarily short 
range and the system has some finite size. Thus, as in 
Rcf. [29|, we introduce a short-range cutoff, e, and a long- 
range cutoff, R. The Fourier transform can then be per- 
formed to obtain 
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where ^2(2;) — {ix^ — l) /2 is the second Legendre poly- 
nomial, 9q is the angle between q and E, and ji [x) is 
the spherical bessel function of the first kind. 

If there is strong confinement along z, then assuming a 
fixed gaussian profile in that direction n {k^) = e^^^"^ 
and integrating Eq. |3] over this momentum profile, the 
effective 2D interaction is 
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II. DIPOLE POTENTIAL IN MOMENTUM 
SPACE 



We consider a gas of one-component fermionic 
molecules with number density n, each possessing dipole 
moment d, aligned by an external electric field E so that 
the inter-molecular interaction potential at large separa- 



where Op is the angle between E and the z-axis, (j)q is the 
azimuthal angle between q and the projection of E onto 
the x—y plane, and we have taken qw — > 0. Alternatively, 
fixing r to lie in the x — y plane and taking the 2D Fourier 
transform of Eq. [T] with a short-range cutoff, w, as in 
Ref. [29] yields the same result with a numerical factor 
of order unity in front of w. 

Similarly, when there is strong confinement and a fixed 
gaussian profile n(fcj^) — e~^^^ in the radial direction, 
integrating out the radial degree of freedom in Eq. [3] 
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yields an effective ID interaction 
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where 9e again the angle between E and the ^-axis, 
T{a,x) is the incomplete gamma function, and 7 « 0.577 
is Euler's constant. (Here again, performing a ID Fourier 
transform of a 1/ potential with a short-distance cutoff 
w gives a similar dependence on q and w, although the 
numerical factors are different.) Unique to the quasi- 
ID case, the momentum dependent part of the effective 
interaction is not independent of the transverse width, 
even for very strong confinement. 



III. UNIFORM SELF-ENERGY, 
COMPRESSIBILITY, AND COLLECTIVE 
MODES AT FINITE TEMPERATURE 



of the dipolar interactions. In the uniform case, to first 
order, n = /uq + S**" (fcFo)j where /zq is the noninteracting 
chemical potential and E"'" (^fo) is the isotropic part of 
the Hartree-Fock self-energy at the unperturbed Fermi 
surface, 
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with no(k) = 1/ \^e^f^^k^/^"^'t^o)/kBi ^ ^j^g noninter- 
acting Fermi distribution. The zero-temperature self- 
energy has recently been calculated in Ref. [1^. Below 
we extend this to finite temperature and obtain the in- 
verse compressibility, effective mass, and zero sound dis- 
persion. 
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1. Compressibility and effective mass 



To calculate the compressibility, n = [3, we 

first need to obtain the chemical potential in the presence 



Let us begin by considering the 3D case at finite tem- 
perature t = ksT /Ep, with /Jo = fio/Ep ~ 1 — 7r'^t'^/l2, 
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FIG. 1; Effective mass at the Fermi surface for KRb in 3D (a) vs temperature at fixed density and (b) vs density at fixed 
temperature, (c) Regions of positive and negative radial effective mass at the poles of the Fermi surface. In all panels, the 
dotted lines show the low T /Tp approximations. 



where 6'k is the angle between k and E, LIq {x) is the 
polylogarithm function, and C' {x) is the derivative of the 
Riemann zeta function. 

In this case, the self-energy has a purely d-wave angular 
dependence and does not affect the compressibility, which 
is given in units of the zero-temperature, noninteracting 
compressibility kq = by 
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The temperature-dependent effective masses along the 
radial and angular directions can be written to first-order 
as in the zero-temperature case (29| : 
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For general temperatures we must evaluate Eq. ([SJ and 
its derivative numerically. At low temperatures, though, 
Eqs. (fTT|) and (fT3|) give analytic approximations on the 
Fermi surface. We have plotted the numerical inverse 
effective mass results at the Fermi surface as functions 
of temperature and density in Fig. [1] as well as the low- 
temperature expansions. 

Note that the deviation of the inverse effective ra- 
dial mass from the inverse bare mass is proportional to 
— P2(cos6'k) (the integrand in Eq. ([5]) is always negative), 
so the effective radial mass is smaller than the bare mass 
at the equator and larger (or negative) at the poles. Also, 
the effective angular mass is proportional to esc 29u- The 
average effective mass over the Fermi sphere is thus un- 
changed by the interaction, but when the effective radial 
mass at the poles of the Fermi surface is negative, the sys- 
tem is unstable. In the low-temperature limit we have 
the stability condition A30 < I/tt, in excellent agree- 
ment with previous numerical results using a variational 
approach which gave X^d < 0.32 27|. Extending this 



result to finite temperature, we obtain the critical line 
for KRb in the n — T plane plotted in Fig. [TJc) . How- 
ever, the system actually becomes unstable at densities 
an order of magnitude lower due to formation of density 
waves perpendicular to the electric field [1^ . 



2. Zero sound mode 

We now turn our attention to the zero sound 
mode, which has recently been calculated in the zero- 
temperature limit [2^ . We anticipate that zero sound 
will only propagate at low temperatures, so we assume 
that the Landau Fermi liquid description 
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remains valid, with the temperature dependence enter- 
ing through the single-particle energy, e(k) = h'^k^/2m + 
S3D (k), and the Fermi distribution function, n(k) = 
1/ (^e'''^^^-t^o)/kBT _|_ r^Yie temperature dependence of 
the Hartree-Fock quasi-particle interaction, /(k, k') ~ 
^30 (q — 0) — VsD (k — k'), is negligible since the inter- 
action near the Fermi surface is independent of the rnag- 
nitude of the momenta in the limit |k — k'| <^ kp [47j . 
Note that, although we have not made the notation ex- 
plicit, the interaction depends on the direction of mo- 
mentum transfer. 

Neglecting collisions, the linearized Boltzmann equa- 
tion [331 is 



6iy (k) = 



q-Vke(k) dn{k) f d^k 



q • Vke (k) de (k) J (27r) 
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where Sv (k) is the deviation of the quasi-particle dis- 
tribution from the equilibrium anisotropic distribution, 
n (k). Integrating out the radial degree of freedom and 
decomposing into spherical harmonics, Yim (f^k)' with 
the angle in the argument measured relative to the direc- 
tion of momentum transfer, q, we can rewrite the above 
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FIG. 2: (Color online.) Zero sound speed (a) and damping (b) vs angle of propagation for finite temperatures. Dashed lines 
show the boundary of the single-particle excitation continuum. 



as 
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with w_Fo = h?kpo/m the noninteracting Fermi velocity and 
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In the preceding we have used the Hartree-Fock single- 
particle energy in order to account for the anisotropy and 
obtain results consistent with Refs. [1^ H^. Since the 
self-energy also appears in the Fermi distribution, the in- 
tegrand is evaluated in some thermally broadened range 
of momenta around the distorted Fermi surface. How- 
ever, note that in order to consistently include only first- 
order corrections to the noninteracting single-particle en- 
ergy, the self-energy terms themselves must be evalu- 
ated around the unperturbed Fermi surface. Hence the 
added superscript denoting the subtle distinction, and 
for k not too far below the Fermi surface one can simply 



take Eg^ (k) = J^^jj (k — kp -I- kpo), where the latter is 
given by Eq. ^ and the three vectors in the argtmient 
are coUinear. 



We will keep only the coupling between the s-wave and 
longitudinal p-wave modes as in Ref. [29|, which yields 
good agreement with the results of Ref. [4^ . The Landau 
parameters for the dipolar interaction have previously 
been calculated (29l . Issj , and in our notation and choice 
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of basis the relevant nonzero parameters are 

/oo,o = 2^714 (cos ^k) , /ll,0 = ^25^4A3DP2(cOS0k) . 

(24) 

The energy, fl, and damping, 7, of the collective modes 
for a given momentum, q, are then given by the solutions 
of det \I - M (q, - i7)| = 0, where 

M (a ui\ — \ '^'^OO'™ -^OO'" Xio,oo (q, i^) /ii,o \ 

V^io.oo (q, /oo,o xii.oo (q, /ii,oy ' 

(25) 

We numerically find solutions corresponding to under- 
damped zero sound propagation. In that case, the dis- 
persion is given by 

Redet|/-Af(q^O,f7^uo(6'g)'7)| =0 (26) 

and the damping, assumed to be small compared to f2, 
is given by 

^ Imdet|J-M(q^O,r!^t;o (gg)g)| 
^ £L=„„(.,),Redet|/-A/(q^O,c.)| 

The dispersion and damping for various temperatures are 
shown in Figs.[2ja) andllfb), respectively, with the den- 
sity chosen such that \zd — I/tt^ for comparison with 
Refs. m, HI]. For temperatures up to about O.lTi?, the 
thermal effect is to increase the axial propagation speed 



while decreasing the range of angles for which zero sound 
propagates. As the temperature increases further, the ax- 
ial propagation speed decreases. Note that at finite tem- 
perature, the mode can become overdamped long before 
it enters the single-particle excitation (SPE) continuum, 
and Eq. ([26|) ceases to have a solution. For temperatures 
above 0.2Tf the damping is large even at 0q = 0, and 
above ~ 0.23rF there is no solution to Eq. (|26p for any 
Qq. Clearly, the assumption of 7 ^ quickly breaks 
down with increasing temperature and propagation off 
the z-axis. 



B. 2D 

1. Compressibility and effective mass 

In the remainder of this manuscript our primary focus 
will be on dipoles aligned perpendicular to their motional 
degrees of freedom, so it will be unnecessary to distin- 
guish between distorted and undistorted Fermi surfaces, 
and Icfo, and we will drop the extra subscript. 

For the 2D case, with fio ~ t\n (e^/* ~ l), 6e the angle 
between E and the normal to the plane, and the an- 
gle between the Fermi wavevector and the projection of E 
onto the plane, from the isotropic part of the interaction 
we obtain: 
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(30) 
(31) 

(32) 
(33) 

m 

m 



with K (x) and E {x) the complete elliptic integrals of the first and second kind 'S] . Although we will focus on 
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the most stable scenario of dipoles aligned perpendicular the self-energy arising from the anisotropic part of the 
to the plane, for completeness we also give the form of interaction: 



E^S (k) = -2^\2D sin2 / -^no(fc') / d0fe. 
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In Eqs. (150]) and ([57)) we have noted that at low temper- 
ature the integrand in the temperature-dependent term 
is only appreciable near x = 1 and expanded the elliptic 
integrals about that point, permitting analytic integra- 
tion. We have shown the derivative of the self-energy 
with respect to density since the inverse compressibility 
For arbitrary temperatures, the com- 



-1 - ^Z^ii 

dn 



is K 

pressibility is obtained via numerical integration from 
Eq. For low temperatures, recognizing that the 

Fermi energy, Ep, and reduced temperature, i, and cou- 
pling strength. A, all depend on density, we can obtain 
an analytic result by differentiating the self-energy given 
in Eq. ([5T|) according to 



K 



dEp 



d 



t d 



A d 



dEp Ep dt 2Ep dX 



'4 vr^ 

= 1 + 16v^A2dP2(cos0b) ( 3 - 



(kp) 

(39) 

+ 0{t'^\nt) . 

(40) 



When the orienting electric field is perpendicular to 
the plane {Op = 0), the compressibility is always posi- 
tive (assuming the gas is not so dense that the ap- 
proximation to the potential breaks down) . As shown in 
Figs.[3Ja) andini^b), the compressibility is nonmonotonic 
in temperature for experimentally relevant densities. The 
nonmonotonicity is due to the interaction-dependent 
term, [kp] /dn, which decreases quadratically with 

t before turning over and increasing with t for fixed den- 



sity. Since t ~ mT/n and A2D 



both the 



temperature at which the peak occurs and the height 
of the peak are larger for higher density. In Fig. ^c), 
we show how the compressibility varies as the density is 
changed. Note that we do not use the natural unit of 



Ko = 2^^"n^ this plot, since it is density-dependent. 
In the low-density (high-i) limit, the compressibility be- 
haves classically, as expected, and is drastically reduced 
in the high-density (low-t) limit. 

The nonmonotonicity in the compressibility shown in 
Figs. |31[a) and |31[b) occurs well below typical experi- 
mental temperatures of ~ 200 nK. However, we may 
also consider the ratio of the finite-temperature interact- 
ing and zero-temperature noninteracting compressibili- 
ties, k/kq ~ n^K ~ dn/d^, which has the same tempera- 
ture dependence as the compressibility. Figures [Sjc) and 
show that this quantity also displays a nonmono- 
tonic density dependence. Moreover, this nonmonotonic 
behavior is evident even at relatively high temperature, 
with T ^ Tp. [For the homogeneous system considered 
here, Tp = 24nK (240nK) for n = 10^ cm'^ (10^ cm'^).] 
In this scenario, the peak is not due to nonmonotonicity 
in the interaction-dependent term, dY,^^° {kp) /dn. Both 
{kp) /dn and d^o/dn are monotonic in the density, 
but while dE"'" (^f) /dn increases with density, dfio/dn 
decreases. The competition between the interacting and 
noninteracting density dependencies is the source of the 
nonmonotonicity. 

Figure [3lje) gives a sense of how k/kq changes as tem- 
perature and density are simultaneously changed in some 
manner. Although, it is hard to see on the logarithmic 
scale, the vertical portions of the lines actually bulge out- 
wards slightly, consistent with the fixed density results 
plotted in Figs. [2Ka) and[3fb). Clearly, nonmonotonic 
behavior is more pronounced if the measurements are 
performed at fixed temperature and variable density than 
the other way around. We will discuss measurement in 
Sec. El 

With the radial and angular effective masses, m* and 
TO^, defined as in Eqs. ([TS]) and pB)) . we have plotted the 
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FIG. 3: Compressibility for KRb in 2D for 6e ~ 0, corresponding to dipoles aligned perpendicular to the plane, (a-b) k/kq vs 
temperature at fixed density; (c-d) k/ko {i.e., dn/djj.) vs density at fixed temperature; (e) k vs density at fixed temperature; 
(f) contour plot of k/kq vs temperature and density (lighter is higher). 



radial mass at the Fermi surface for the isotropic case in 
FigHl We see that at low temperatures or high densities 
the effective mass is diminished by the interaction while 
at high temperatures or low densities the effect of the 
interaction is small and the effective mass approaches its 
bare value. The analytic low temperature results from 
Eq. ((55)) are shown as dotted lines. 



2. Zero sound mode 

For dipoles aligned perpendicular to the plane, the 2D 
interaction is isotropic in contrast to the 3D case, and 
we need not include coupling to higher partial waves 
when calculating the zero sound. Likewise, since the 
single-particle energy is also isotropic, to leading order 



we can simply use the bare single-particle energy. The 
zero sound dispersion is then simply given by the solution 
of 1 — V2D (0) X (9 ^ 0, w voq) = 0, where 



d^k dn{k) q- 



X (g ^ 0, Voq) = 

xdx — e( 



:>)/t 



27r 



t [e(^'-Ai'o)/* + 1] 



1 - 



VQ 
Vp 



(41) 



Note from Eq. that the thickness of the monolayer 
plays a critical role in determining the zero sound speed. 
The dispersion is plotted in Fig. [5] as a function of tem- 
perature and density for a thickness of w = lOnm. 
In the limit T/Tp — >■ 0, one can show that vo/vp = 
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FIG. 4; Effective mass at the Fermi surface for KRb in 2D 
(a) vs temperature at fixed density and (b) vs density at fixed 
temperature. Dotted lines show the low T /Tf approxima- 
tions. 



y/2dhnJi^JwwW |49i], and we numerically recover this 
limit to better than 2%. The damping is negligible, and 
we have not shown it. The mode propagates even at 
quite high temperatures due to the strong repulsive delta- 
function core of the effective 2D interaction (U). In fact, it 
persists to arbitrarily high temperature for tight enough 
confinement. This is an artifact of the ultraviolet di- 
vergence of the 2D Fourier transform of the interac- 
tion, and in experiment this robustness will be limited by 
the range, r*, at which the actual intermolecular poten- 
tial deviates from a behavior. However, provided 
kprsf <C 1 and <C w, these results will remain valid 



and the collective mode will be robust at experimentally 
feasible temperature and density. 

n=10^ cm~^ 




50 100 150 200 

T(nK) 
T=200nK 




n (cm ^) 

FIG. 5: Zero sound speed for a f nm thick monolayer of KRb 
(a) vs temperature at fixed density and (b) vs density at fixed 
temperature. 



C. ID 

1. Compressibility and effective mass 

For the ID case, /Iq ~ 1 + 7r^/:^/12, and 9e is the 
angle between E and the trap axis. In this case, recall 
from Eq. ([5]) that the results will depend explicitly on 
the transverse width (or short range interaction cutoff), 
w. In the limit kpw"^ ^ 1, the Hartree-Fock self-energy 
is 
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FIG. 6: Compressibility for KRb at 6e = 7r/2, corresponding to dipoles aligned perpendicular to the axis in quasi- ID. (a,d) 
k/ko (i.e., dEp/dy?) vs temperature at fixed density; (b,e) k vs density at fixed temperature; (c,f) k/kq vs density at fixed 
temperature; (g,h) contour plot of k/kq vs temperature and density (lighter regions are higher). 



SiD(k) = — Aiz5-fP2(cos0B) / dk'no{k'){k-k'f {-f + 2\n\{k-k')w\) 



(42) 



SiD (kp) = —\iDEFP2{cOSeE) 



— XidEfP2{cos Oe) 



(7 + 21n(fcFw)) I 2 - ^f^' uif-e-^'"/* 



- ( 7 - - + 2 In {2kFw) ) + (7 + 1 + 2 In {2kFw)) —t 



dx (1 — x) 



(9 [t^ Int] 



Inll -xl 



(43) 
(44) 



Differentiating with respect to density yields 

8(7 + 21n(2fci.w)) 



9 2 

— = l--77t' + ^AiDP2(cOSf?B) 



(7-l + 21n(2A:j.w))^t2 





0[t^\nt\. (45) 



However, even if we take transverse size w ^ lOOnm, 
corresponding to a very strong optical dipole trap, then 
at n ^ 10^ cm^^ although the system is quasi- ID 
{h?/mw'^ > kBT,EF) for T < 40 nK, fc^w is actually 



of order unity. This tends to diminish interaction effects 
since less repulsive off-axis interactions must be included 
in the effective ID interaction. (We did not hesitate to 
take kFW <C 1 in the 2D case @ because there the work- 
ing assumption is that several independent quasi-2D lay- 
ers will be formed by an optical lattice, which can be 
made considerably tighter than a dipole trap. In any 
event, the general treatment of the multilayer quasi-2D 
in Sec. II VI includes finite thickness effects.) For arbitrary 
values oi kFW, then, a similar calculation without assum- 
ing qw m Eq. ([5]) yields 
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0{t^lnt) (46) 



dEiD (fc) 



dk 



k—kp 



12 



(1 + 20fc>2 + 324u;4) 6^'=^'" r(0, 4kj, 



7 - (2) + 84u;2 - In kj^w^ 



^t^\nt + 0{t*lnt)^ (47) 



Ko 



' ^XiDP2{coseE) 1 86^'=-'" r(o,4fc>2) 



6 



+8klw^ { 1 + e'"=J-'"'r(-l, 44^2^ ) +6G'2;3 | 4^ 



1 _ e4fc^™ p^Q^ 44^2) (3 + 4fc|w2 + 324 
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-3/2,-1 
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r 



-O {t^ Int) 

(48) 



where G^'^ ^ 



2,3 y-^ bTb2b3) ^ Meijer G-function [48|. The 
more general Eq. (j48l) was used for the low temperature 
approximation in Fig. |6l For the second row of plots, 
with w = 10 nm corresponding to a tight optical lattice 
forming an array of independent quasi- ID tubes, the sim- 
pler initial approach gives the same results except for at 
high densities. 

Figures |6l[a) and |6jd) show a maximum in the com- 
pressibility at finite temperature as in the 2D case. Now, 
though, t ^ mT/n^, so the temperature at which the 
peak occurs increases with density more rapidly than in 
2D. For w = lOnm, we have plotted results for higher 
densities and temperatures since these values can be 
reached while maintaining low-dimensionality. In con- 
trast to the 2D case, in a tightly confined system at rea- 
sonable density, the peak occurs at higher temperature as 
shown in Fig.lH^d). However, it is important to note that 
in ID this is no longer a signature of interaction effects 
since the compressibility exhibits a maximum even in a 
noninteracting system due to the nonmonotonic temper- 
ature dependence of dfio/dn. In fact, the main effect of 
interactions is to add a constant shift to the inverse com- 
pressibility, making the peak relatively less dramatic. 

As in the 2D case, we can also look for nonmonotonic- 
ity in the density dependence of k/kq — dEp /dfi ^ tv'k. 
In the 2D case, we found a maximum in k/kq in the pres- 
ence of interactions. Here there is a weak maximum in 
the absence of interactions, due to the nonmonotonicity 
of In the presence of interactions, this feature washes 
out, but the interactions may give rise to a more pro- 
nounced local maximum. Unlike the 2D case, here the 
nonmonotonicity in density is predominantly due to non- 



monotonic behavior of c?S {kp) j dn rather than competi- 
tion between the interacting and noninteracting terms. 
For dipoles perpendicular to the trap axis and w ~ 100 
nm, the density dependence at fixed temperature above a 
few uK is monotonic as in Fig. IHl^c). However, for tighter 
radial confinement, the interaction effects become more 
pronounced and nonmonotonicity emerges at higher tem- 
peratures, as shown in Fig. [6l^f). At high densities, the 
effective ID interaction ([5]) at the Fermi surface saturates 
due to the dependence on the finite transverse size. As a 
result k/ko begins rising again. Actually, for w — lOOnm 
(w = lOnm), one-dimensionality is already beginning to 
break down for n > 2 x 10* cm^^ (n > 2 x 10^ cm^^), 
since there Ep ^ jmrn^ . 



The effective mass is plotted in Fig. [7] as a function of 
density and temperature. The effective mass is again di- 
minished by the interaction and slowly increases towards 
its bare value at high temperatures. At low densities, as 
in 2D, the reduction of the effective mass by the inter- 
action increases with density. However, as the density 
increases further the effective mass eventually returns to 
its bare value since the effective ID interaction saturates 
while the Fermi velocity continues to increase. In order 
to clearly show the trend, we have included in Fig. [TJb) 
extremely high densities. The analytic low temperature 
results from Eq. (IT71) are shown as dotted lines. 
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FIG. 7: Effective mass at the Fermi surface for quasi-lD KRb 
with w = lOnm (a) vs temperature at fixed density and (b) 
vs density at fixed temperature. Dotted lines show the low 
T/Tf approximations. 



FIG. 8: Real (solid) and imaginary (dashed) parts of the zero 
sound speed for quasi-lD KRb with transverse size lOnm (a) 
vs temperature (note the log scale) at fixed density and (b) 
vs density at fixed temperature. 



2. Zero sound mode 



As in the 2D case, we consider the zero sound disper- 
sion given by 1 — Vm (0) x — > 0, a; — voq) — 0, where 



Xiq-^ O,cog) 



kp 



dx 



2TTt [g(l.2-Mo)/* + 1]^ 
1 



vf 



(49) 



Once again, the finite transverse size which appears 
in Eq. ([5|) is important in determining the zero sound 
speed. The dispersion is plotted in Fig. [8] as a func- 
tion of temperature and density for w = lOnm and 
dipoles aligned perpendicular to the axis of the tube. In 
the low-temperature limit, one can show that vq/vf = 
^ (P-mj^'KuP'W-kp^ and we numerically recover this limit 
to better than 1%. One cannot obtain this limit by going 



to high density instead of low temperature, since, unlike 
in 2D where x approaches unity in the high density limit, 
here x goes to zero like 1/n for vq ^ vp. Thus, at high 
density the zero sound speed decreases to unity. The 
damping strength, 7/9, is shown as a dashed line and it 
is negligible except near where the mode vanishes. The 
collective mode is even more robust than in the 2D case 
since the ID Fourier transform of the interaction is even 
more divergent. Again, these results will remain physical 
until becomes comparable to either of the other length 
scales, kZ^ or w. 
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FIG. 9: (Color online.) Compressibility for KRb in a periodic potential along z. (a-b,d-e) k/kq vs temperature at fixed density; 
(c,f) k/kq (i.e., dn/dfj,) vs density at fixed temperature. 



IV. 3D MULTILAYER SELF-ENERGY AND 
COMPRESSIBILITY AT FINITE 
TEMPERATURE 

A. Compressibility 

Now we consider a 3D gas placed in a strong ID op- 
tical lattice of depth Vq and wavelength A, resulting in 
a "stack of pancakes" configuration. In the limit where 
Vq and A become infinite, this is the same as the strictly 
2D monolayer considered in Sec. IIII Bl Generally, how- 
ever, we find that the nonmonotonic behavior discussed 
in Sec. IIII Bl persists in the presence of multilayer and 
finite thickness effects. We will restrict our discussion 
to the case of dipoles aligned along the direction of the 
optical lattice. 

Because of the external periodic potential along z, the 
effect of the interaction on the chemical potential to first 
order, fi = iiq + A/i, must now be written in real space: 



= y" d^xd^xV*^ (x) S (x, x') (x') ^ 



(50) 



where the angular averaging over the Fermi surface de- 
noted by the brackets is necessary to separate the shift 



in chemical potential from the anisotropic deformation of 
the Fermi surface, and 



S (x, x') = / d^^" f -^no (k) V3D (x - X' 
J J (Stt) 



: (x")r S (x - x') - 0* (x') (x) 6 (x" - x') 



(51) 



The unperturbed single-particle wavefunction can be 
written in terms of the reciprocal lattice points, which 
are integer multiples of 27r/A, in the Bloch form 

N 

0k (x) = e*-^ "9 (^-) e-*2-9-./A, (52) 

q=-N 

The coefficients Uq (k) are easily obtained numerically, as 
are the corresponding single-particle energies needed to 
evaluate no (k) . In the above we have introduced a cutoff, 
2Nt:/X, on the Bloch momentum in order to perform the 
sum numerically. After some algebra, Eq. (|50p can be 
recast as a sum over momenta. 
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A/i = 



o no (k) > 



0,— j ^ Uq, {k^)u*q,^q{k^) ^ V (fcj^J Wg-'+g (^J^J 
, g"=g< 



g =9< 



g> 



^ Mg' (fcz) Mg'+g (fcfj 



g =g< 



(53) 



where Vsd ((ZijQz) is the 3D Fourier transformed inter- 
action defined in Eq. ^ and g< = max(— A^, —N — q), 
g> = min {N,N — q). We typically take N ~ 50, which 
we find to be suSicient for the lattice depths treated. 

To facilitate comparison with the strictly 2D case, we 
measure density in terms of the 2D density per layer, 
"2D — n^o^, and compressibility in units of kq = 
where E2D = fi.^47rn2i:i/2m. By first analyt- 



1 dn^D 
dE2 



ically performing the azimuthal integration in Eq. (j53p 
and differentiating with respect to density, then numer- 
ically summing over reciprocal lattice points, averaging 
over the Fermi surface, and performing the remaining 
two-dimensional integral, the compressibility can be effi- 
ciently obtained. In Fig. [S] we show the compressibility 
for KRb in a deep optical lattice (here the recoil energy 
Eji = {2tt / Xf / 2m) . For strong lattices, the com- 
pressibility behaves similarly versus temperature to the 
2D results of Fig. [31 with a maximum at nearly the same 
temperature. For weak lattices, the effect is less pro- 
nounced. 

As density is varied, the ratio k/kq behaves as in the 
2D case for low densities, with nonmonotonicity evident 
even for T Tp, but at high density E2D becomes larger 
than the lattice bandgap and quasi-two-dimensionality 
breaks down. Here the compressibility exhibits discontin- 
uous jumps as the chemical potential enters higher bands, 
reminiscent of the effect seen in Ref . [53| • At very large 
densities, the recoil energy is negligible compared to the 
kinetic energy of the gas, and k/kq should increase with 
density as in the 3D case rather than decreasing with 
density as in the 2D case [s^. So, in addition to the 
interaction-induced local maximum analogous to the 2D 
case, there can also be a local minimum as the system 
crosses over to three-dimensional behavior. This is sim- 
ilar to the finite confinement effect seen in the ID case 
above. 



B. Collective modes 

For the 3D multilayer system described above we now 
calculate the low-lying collective modes at finite temper- 
ature, neglecting coupling between partial waves. In the 
limit of a large bandgap, such couplings should be sup- 
pressed, and outside this limit the description should still 
be qualitatively valid. However, because of the spatial 
inhomogeneity, the 3D polarizability can no longer be 
written as a scalar function of the momentum transfer. 
In terms of the Bloch wavefunctions defined in Eq. (|52p , 
the interaction and the zeroth order polarizability can be 
written in momentum space as 



Vk,k' - J d^^£^'<i>l (x) y (x - x') 0k' (x') 



/ X ki,k2* ki,k2 '^o (ki) — no (k2) 

Xk,k' (w) = 2^ a^^' ^ a. ^' ^ 



(54) 



ki,k2 



c^-eo(ki)-f eo(k2)+iO' 
(55) 

where cq (k) is the bare single-particle energy associated 
with </)k, and 

al^ '^- = j £xd^^'<j>l^ (x) </.k2 (x) (/-k (x) . (56) 

The collective modes are given by the zeros of the 
determinant det|l — M(w)| = 0, where Mk^k' (w) = 
{Vx ('^))k k'- Conservation of momentum in the x — y 
plane and of quasi-momentum along z reduces M to a 
matrix function of the momentum transfer, Mm,n (q, i^), 
where m is an integer denoting a reciprocal lattice vector 
Km — 27rni/A and q^ is in the first Brillouin zone. Some 
algebra results in a computationally feasible expression 
for the matrix elements. 
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FIG. 10: (Color online.) Zero sound speed vs density for KRb in a periodic potential along the z-axis. 



3N-m 



Mm,n (q, W) = ^ / 



d^k no (fc_L - q±,kz - + 2ttQ/\) - rip {k±, kz) 
{2nf uj-e{k^- q±, k, - q, + 2t:Q / X) + e (/cj^, fc^) + iO 



XI ( - 9z + ^ ) Ura' (k,) <, ( q^ + 



27rm 



A 

/, 27rg\ , , / 27rn\ ^ / 27r(n-n')' 

2^ WQ+«+rn'-n' ( - 9^ + 1 (Kz)u„' I + — I I (J J. , (J^ + 



(57) 



The matrix elements can be evaluated numerically for a 
given frequency and momentum. To proceed, one must 
obviously truncate the matrix at some finite reciprocal 
lattice vector, and we have used a 17 x 17 matrix, as this 
seems to include the dominant terms for large potential 
depth, Vq- As Vq is decreased, the off-diagonal terms 
of the matrix vanish while the number of nonnegligible 
diagonal elements increases. 

In Fig. [TU] we plot the lowest-lying underdamped col- 
lective modes for long-wavelength, in-plane momentum 
transfer, (?2 = 0, q_L — >■ 0. The zero sound speed is shown 
in units of the 2D Fermi speed, vp = h^^/Aww2D/m. In 
addition to the Landau damping due to single particle 
excitations in the lowest band for vp < vp, there are also 
regions of damping due to single-particle excitations in 
excited bands which occur for vq > vp- The boundaries 
of these regions are shown in dotted lines. The modes 
plotted are weakly thermally damped. 

It is instructive to first consider the 20nK curve in 
Fig. Iior a). where one can identify three distinct regions. 
At densities below ~ lO^cm^^, the thermal energy is 
comparable to the 2D energy scale, E2D — fi^4:TTn2D /2m. 
For higher densities, one might expect to see a density 
independent speed, as in Sec. lIII B 2l and the zero temper- 
ature calculation of Ref . 49] . This appears to be the case 



for an intermediate region, with the value wq > wf ~ 1-4 
in qualitative agreement with the value ~ 1.6 obtained 
using the multilayer result of Ref. [i^ and inserting an 
estimate for the transverse size in the effective 2D interac- 
tion (m based on the harmonic approximation to the local 
wells of the optical lattice. However, for even higher den- 
sities the speed abruptly drops towards the continuum. 
This begins happening where E2D becomes comparable 
to the bandgap, signalling a transition to 3D behavior. 
As we have seen in Sec. lIII A 21 the mode with momentum 
transfer perpendicular to the electric field must vanish in 
3D. In this region, the behavior is roughly independent 
of T, since the reduced temperature T/Tp is effectively 
zero. 

For T = 200nK, the mode does not propagate except 
at high densities, and we cannot reach the low temper- 
ature limit without increasing the density into the 3D 
transition. For the weaker lattice in Fig. [TOlb). even at 
T — 20nK the intermediate plateau has vanished because 
higher bands are already becoming active as we increase 
the density out of the large T/Tp regime. In Fig. ITUTb) 
we see multiple modes entering the continuum, although 
only one is stable at low density. Also, at high density, 
faster zero sound modes appear above the first excited 
SPE band. 
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FIG. 11: Sloshing frequency vs density for KRb in a periodic 
potential along the z-axis. 

In Fig. [TT] we show the lowest underdamped collec- 
tive excitation for momentum transfer along the z-axis 
in units of the frequency of the local harmonic well ap- 
proximation to the optical lattice. In a harmonic trap 
this is the Kohn (or sloshing) mode jsi'l, and it is inde- 
pendent of the interaction and the momentum. Here the 
density dependence is due to anharmonicity of the poten- 
tial. The curve shown is independent of temperature and 
momentum. For Vq = 2bEji, we do not find this mode. 

So, although one can observe a Kohn-like collective 
mode along the direction of the optical lattice at cur- 
rently attainable temperatures and densities, the obser- 
vation of zero sound in a reasonably strong optical lattice 
with standard site separation on the order of a micron 
requires an order of magnitude improvement in the at- 
tainable temperature or density. Based on the 2D calcu- 
lations of Sec. IIIIB 21 it is reasonable to expect this re- 
quirement could be eased by using a more tightly spaced 
lattice. 



V. EXPERIMENTAL DETECTION OF 
COMPRESSIBILITY 

We have shown that a quasi-low-dimensional exper- 
iment may be able to observe an interaction-induced 
maximum in the ratio of the finite-temperature inter- 
acting compressibility and the zero-temperature nonin- 
teracting compressibility, particularly if the measure- 
ment is performed with variable density. Now, the re- 
sults we have shown are for homogeneous systems, while 
there is always a global confining potential in experi- 
ments. However, if we incorporate a weak harmonic trap 
through a simple local density approximation (LDA), 
/i (r) « /i (0) — TOCj^r^/2, where one can think of a local 

compressibility, K (r) n-^ (r) dfi{r) n'^ {r)'muj'^r dr ' 

then the trap can be beneficial in that precise in-situ 



measurement of the density profile automatically allows 
the simultaneous measurement of compressibility across 
a wide range of densities. This technique has recently 
been used to measure local compressibility in an optical 
lattice jS^]. The ratio k (r) /ko (r) ^ dk'jp (r) /d^x (r) is 
similarly accessible. In Fig. [12] we have shown the ratio 
as a function of radial distance in the LDA, and the asso- 
ciated density profile, for KRb at reasonable temperature 
and peak density. A gaussian profile with the same peak 
density and number of molecules as the calculated profile 
is shown in Figs. [T^b) and [T^d) for reference. The 2D 
results shown are generated from the strictly 2D calcu- 
lations of Sec. nil Bl rather than the more complicated 
multilayer case. 

If the strictly 2D case is realized by applying in- 
tense counter-propagating beams of wavelength ^ l/im 
along the axis of a cigar-shaped cloud that contains 

10^ molecules and is 600/ini in length, with radius 
- lO^m 0,i|, then a typical 2D plane will contain ~ 100 
molecules. As shown in Fig. [T2T a). to reach suitably large 
densities with such a small number requires the in-plane 
trapping frequency to be strong enough that the spatial 
extent of the cloud is rather small. This may preclude 
in-situ measurement of the profile, so for that specific 
method, it is preferable to apply the lattice perpendic- 
ular to the axis of the dipole trap in order to deal with 
larger clouds. Rather than a cigar-shaped dipole trap, 
one might instead apply a pancake-shaped trap as in 
Ref. [ll|, but it would be rather difficult to reach the 
2D regime without using an optical lattice. 

If the ID case is realized by applying two orthogo- 
nal pairs of counter-propagating beams perpendicular to 
the axis of the initial cigar-shaped trap, then a typical 
quasi-lD tube may contain ^ 1000 molecules. This case 
is shown in Figs. [T2jc) and FlW d'). In addition to the 
interaction-induced maximum, there is a local minimum 
for sufficiently high densities, as discussed above. For 
those higher trap frequencies, the density decreases more 
rapidly and the features are correspondingly narrowed. 

However, there is the possible complication that in the 
presence of a trap, the Hartree term is no longer just a 
constant. Due to the long-range interactions, one may 
need to explicitly take the nonuniformity into account 
in a self-consistent manner with a spatially dependent 
Hartree term, as was done, for example, in Refs. |26l. [53| . 
It is not clear how important this effect is, although one 
may estimate that it will be negligible if the interaction 
energy, d'^/£f, at distances on the order of the character- 
istic length scale of the trap, £t, is very small compared 
to the local Fermi energy. Since we have already assumed 
the local Fermi energy is much larger than h? /2mi1 by 
taking the LDA, this implies 2md'^ /h? ^ is a sufficient 
condition for neglecting spatial variation of the Hartree 
term. For KRb, 2md^ /h? ~ 30nm, so it is reasonable 
to treat the system the way we have as a first step. Fur- 
ther investigation is required, though, to fully resolve the 
issue. 
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FIG. 12: 2D results shown in panels (a,b) and quasi-lD results with w — 10 nm in panels (c,d) for KRb. (a,c) Local k/ko vs 
distance from trap center, using LDA. (b,d) Density profile, using LDA, with a gaussian profile shown for reference. 



VI. CONCLUSIONS 

In this work, our goal has been to propose well-defined 
experiments which can see quantum many-body effects in 
a stable, ultracold gas of dipolar molecules at finite tem- 
perature. Our work is admittedly a leading-order approx- 
imate calculation, however it should be generally valid 
qualitatively and, for the coupling strengths accessible to 
experiment, essentially exact. Our results indicate that 
the compressibility of a quasi-low-dimensional gas should 
exhibit a nonmonotonic temperature dependence for low 
temperatures. (Note that a similar maximum in the 
Hartree-Fock compressibility occurs in an electron gas 
with 1/r interactions, although there it is due to at^lnt 
term in the self-energy [sj]; see Appendix A.) Further- 
more, the closely related quantity dEp/dii, which is the 
ratio of the finite-temperature interacting compressibility 
and the zero-temperature noninteracting compressibility, 
should also exhibit a quantum many-body effect, namely 
a nonmonotonic density dependence, which persists for 
T ^Tp. 

We have also considered the finite temperature behav- 
ior of the effective mass and zero sound dispersion. In 
the 3D case, the zero sound shows a nonmonotonic tem- 
perature dependence, in addition to angular anisotropy. 
However, as expected, we have seen that the damping 
increases rapidly with temperature, so that the mode 
does not propagate at currently feasible temperatures 



and densities. However, in quasi-low-dimensional scenar- 
ios, the mode is far more robust, propagating at relatively 
high temperatures for tight transverse confinement, as- 
suming the length scale at which the interaction deviates 
from behavior is small compared to the transverse 
size. 

We have performed a numerical calculation for a 3D 
dipolar gas in a ID periodic potential, and in addition to 
finding similar behavior as in the 2D case in the presence 
of multilayer effects, we observe effects of excited bands 
at high density and the breakdown of two-dimensionality. 
In particular, at moderate lattice depths we see discon- 
tinuous jumps in the compressibility and interesting col- 
lective mode structure. 

Finally, we have discussed a method for experimentally 
observing the characteristic nonmonotonicity of the com- 
pressibility in the presence of a global confining potential 
via high-resolution in-situ imaging. 

JFK acknowledges helpful discussions with D.-W. 
Wang. This work is supported by JQI-AFOSR-MURI. 



Appendix A 

Below we briefiy compare the results given above for a 
spinless 2D dipolar gas to results similarly derived for a 
2D spinless Coulomb gas with V (r) = e^/r. In the latter 
case, the low-temperature behavior within Hartree-Fock 
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is Hi], 

^„l + l(^l + 0,lM' + ^t'lrit] (Al) 

K TT \ 32 / 

where = l/as^Trn, with the Bohr radius. The 
In t term is not an artifact of the dimensionahty, as a 



similar result has also been obtained for the 3D Coulomb 
gas [13] • Here again, the compressibility behaves non- 
monotonically, due to the interaction-dependent term 
first decreasing, then increasing with temperature. How- 
ever, in contrast to the dipolar gas, this nonmonotonicity 
is evident already at leading order in t. 
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